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In contrast to a complex feature of antinodal 
state, suffering from competing order(s), the 
“pure” pairing gap of cuprates is obtained in the 
nodal region, which therefore holds the key to the 
superconducting mechanism. One of the biggest 
question is whether the point nodal state as a 
hallmark of d— wave pairing collapses at T c like the 
BCS-type superconductors, or it instead survives 
above T c turning into the preformed pair state. 
A difficulty in this issue comes from the small 
magnitude of the nodal gap, which has been pre¬ 
venting experimentalists from solving it. Here we 
use a laser ARPES capable of ultrahigh energy- 
resolution, and detect the point nodes surviving 
far beyond T c in Bi2212. By tracking the tem¬ 
perature evolution of spectra, we reveal that the 
superconductivity occurs when the pair breaking 
rate is suppressed smaller than the single parti¬ 
cle scattering rate on cooling, which governs the 
value of T c in cuprates. 

In cuprates, the energy gap (pseudogap) starts open¬ 
ing at a temperature much higher than T c , in some cases 
above room temperature. Many experimental evidences 
[1-9] point to a competing-order origin, rather than the 
preformed pair, for the pseudogap observed around the 
antinode with the maximum energy gap. The pseudogap 
state has been revealed to sharply diminishes toward the 
antinode and disappear far off the node by the spectral 
weight analysis for the ARPES data of optimally doped 
samples [9]. This is compatible with the recent results 
by the scanning tunneling spectroscopy (STM) and the 
resonant x-ray scattering, showing that the density wave 
Q-vectors are detected only in the Cu-0 bond directions 
[5, 7, 10-12]. While the whole Fermi surface is eventually 
dominated by the pseudogap close to the Mott insulat¬ 
ing phase [6, 13, 14], the “pure” pairing state seems to 
be realized around the node at least in the optimally- 
and over-doped regime [9, 15-21]. The relevant feature 
to it would be the spatially homogeneous electronic state 
seen at the low-bias in the STM spectra, which reflects 
the nodal momentum region [11, 20, 22]. It strongly con¬ 
trasts to the highly inhomogeneous spectra at the high- 
bias associated with the antinodal states. Unveiling the 
nature of the spectral gap near the node is therefore 


crucial to elucidate the superconducting mechanism in 
cuprates. A difficulty however is the small magnitude of 
the gap, which has been challenging the experimentalists 
to investigate. 

It has been proposed that the pairing-gap evolution 
with temperature simply follows the conventional BCS 
function [15], and Fermi arcs (disconnected segments 
of gapless Fermi surface) [23] emerge at T c [6, 15- 
17, 20, 24], marking momentum borders between the 
superconducting and the competing pseudogap regions 
[6, 15]. However, it seems contradicting the observations 
of the Nernst and diamagnetic effects above T c [25, 26], 
which are viewed as signatures of a phase incoherent su¬ 
perconductivity. 

Recently, a contrasting view was proposed [9, 18, 19]: 
its underlying idea is that one should discard the notion 
of electron quasiparticles, instead pay attention to the 
density of states, which is an effective way of judging 
the existence of energy gap. Accordingly a momentum 
integration of angle-resolved photoemission spectroscopy 
(ARPES) spectra has been performed over a selected part 
of the momentum space. This quantity contributed from 
the nodal region was found to have a gap-like structure 
even above T c [18, 19]. The result seems to be in direct 
opposition to the above widely accepted view. Neverthe¬ 
less, the evidence for single-particle gap with the point 
nodes surviving above T c is still missing, and strongly 
desired to unveil the nature of Fermi arc, which is tiled 
to the pairing mechanism of cuprates. 

The determination of the momentum-resolved gap 
structure has been also attempted by STM studies 
through the spectral line-shape analysis [20], and by 
applying the octet model to the interference pattern 
[21, 27, 28]. While these STM techniques are very suc¬ 
cessful, they seem to be limited to the investigation of 
the antinodal region; the gap structure close to the node 
is not determined [27, 28], or the gapless Fermi arc is 
obtained even below T c [20, 21], which has never been 
reported by the ARPES studies. It is thus crucial to in¬ 
vestigate the nodal region by means of the ARPES, which 
is specialized for an observation of the momentum space. 

Here we examine the momentum-resolved single¬ 
particle spectra of Bi2212 obtained by a laser ARPES 
[29]. The ultra-high energy resolution and bulk- 
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sensitivity achieved by utilizing a low-energy laser source 
(hv= 7eV) enabled us to obtain high quality spectra with 
an extremely sharp line shape. We demonstrate, within 
the quasiparticle picture, the absence of the gapless Fermi 
arc at T c , and an isotropic temperature evolution of point 
nodal pairing state persistent far above T c . We find that 
not only the single-particle scattering rate (r s i ng i e ), but 
also the pair breaking rate (r pa i r ) is required to repro¬ 
duce the ARPES spectra. Furthermore, the magnitude 
of T c is determined by the mutual relation between r sing i e 
and r pa i r in the form of T S i ng i e (X c )— T pa j r (X c ). Impor¬ 
tantly, the momentum-integrated spectra of ARPES and 
STM previously investigated are not capable of separat¬ 
ing these two quantities (r s i ng i e and r pa i r ) (see Supple¬ 
mentary Note 4), thus the present results provide novel 
ingredients essential to formulate the pairing mechanism 
of cuperates. 

In Fig.l, we show typical data obtained inside the 
nodal region where the Fermi arc (bold orange curve in 
the inset of Fig. Id) was previously claimed to appear at 
T c [6, 15, 16, 20, 24]. The ARPES intensity map divided 
by the Fermi function (see Fig.la) shows an energy gap 
and an upper branch of the Bogoliubov dispersion at low 
temperatures, as an indication of the pairing state. We 
extract the spectra at the Fermi momentum (kp) over a 
wide range of temperature in Fig.lb, and plot the peak 
energies (£ pea ks) in Fig.le. In the same panel, we also 
plot £ pea kS of EDCs symmetrized about the Fermi en¬ 
ergy (Fp) to remove the effect of Fermi cut-off [23], and 
confirm a consistency between the two different methods. 
Our high quality data clearly exhibit that the gap is open 
even at T c (=92K) (see green spectra in Fig.lb and lc), 
and the £ pea k(T) disagrees with the BCS gap evolution 
(blue solid curve in Fig.le). Even if assuming a phase 
fluctuation slightly above T c , the BCS-type curve (blue 
dashed curve) still does not fit to the data. 

To pin down the cause of this anomaly, we examine 
the momentum variation of £ pea k(T) for the optimally 
and overdoped samples (OP92K and ODT2K) in Fig. 2a 
and Fig. 2b, respectively. Surprisingly, the gap does not 
close at T c regardless of kp points. The symmetrized 
EDCs around the node for OP92K are plotted in Figs. 
3b (T = 10K) and 3c (T = T c ). We find that the d-wave 
gap with a point node persists at T c (Fig.3d); the Fermi 
arc is absent. While a small uncertainty in gap estima¬ 
tion remains in the close vicinity of the node due to the 
finite spectral width, it is negligible compared with the 
previously reported extensive gapless Fermi arc (orange 
arrows in the inset of Figs. 2a and 2b). The absence of 
Fermi arc is further confirmed in Fermi-function divided 
band dispersions (Fig. 3a) measured at T c along several 
momentum cuts (color lines in the inset of Fig. 3d). The 
loss of spectral weight at Fp due to the gap opening is 
seen in all the maps except for at the node (see Sup¬ 
plementary Fig.6 for more details). Our high resolution 
data also show other inconsistencies with the previous 


expectations [6, 15-17]. First, the length of arc with sin¬ 
gle spectral-peaks (£ pea k = 0) is not linearly extrapolated 
to zero at T=0 against the nodal liquid behavior (Sup¬ 
plementary Fig.7) [16, 17]. Secondly, the temperature 
evolution of such an arc is gradual up to far above T c 
with no indication of momentum borders separating two 
distinct states [6, 15]. 

For a further examination, we normalize each curve of 
^peak(T) to the maximum value at the lowest temper¬ 
ature in the bottom panels of Figs. 2a (OP92K) and 
2b (OD72K). One can confirm that the data are mis¬ 
matched with the conventional BCS curve (red dashed 
curves) even in the close vicinity of the node. More im¬ 
portantly, the £ pea k(T) behavior with a steep drop to zero 
becomes more gradual with getting away from the node, 
and it eventually follows a BCS-type gap function (green 
curves) with an onset much higher than T c (~ 135K and 
~ 89K for OP92K and OD72K, respectively). 

To clarify the anomalous feature above T c , here we 
investigate another spectral measure, so-called leading- 
edge of EDC, which is also commonly used for a gap es¬ 
timation. Figure 4 a-d show the non-symmetrized EDCs 
of OP92K, normalized to the peak intensity of each curve. 
The energy location of spectral leading-edge (£le), a t 
which the spectral intensity becomes half (white allows 
in a-d), is plotted in Fig. 4e-h as a function of tem¬ 
perature. The nodal values (Fig.4e) have a T-linear be¬ 
havior, which is expected for the spectra dominated by 
the Fermi cut-off effect. While the nodal £le (T) could 
have a complex behavior [30] with a poor energy resolu¬ 
tion, such a effect seems not to be observed in our case. 
Hence, even the slightest energy gap could be detected 
as the deviation of £le(T) from the T-linear behavior. 
Such a deviation is indeed observed for all the k points 
off the node (Fig.4f-h). It is more clearly demonstrated 
in the bottom panels of Fig.4f-h, where the difference of 
£le(T) from the T-linear behavior is extracted. An as¬ 
tonishing result we found is that the onset temperatures 
are almost the same (~ 135K) regardless of the directions 
(0s), which validates that the point nodal state persists 
up to ~ 135K. 

Importantly, the gap opening temperature is coinci¬ 
dent with the onset of the BCS-type gap function ob¬ 
tained with the symmetrized EDCs far off the node (a 
green curve in Fig.2a). A plausible explanation for it is 
that the spectral peak energy underestimates the “real” 
energy gap (e pea k < A) [9, 31, 32], and energy gaps (As) 
comparable to or smaller than the peak width cannot be 
detected by tracing the peak positions (see a simulation 
in the Supplementary Fig.8). This situation could occur 
at high temperatures, and it gets more serious toward the 
node with smaller A. We actually detect the signature 
of such a small gap in the off-nodal spectra with a sin¬ 
gle peak above T c (Supplementary Fig. 10); the spectral 
width becomes smaller with increasing temperature up 
to ^ 135K, which contrasts to the monotonic broadening 
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seen in the nodal direction. The characteristic momen¬ 
tum variation of £ pea k (T) in Figs. 2 a and 2 b, therefore, 
could be a natural consequence of A (T) having the same 
onset temperature at ~ 135K and ~ 89K, respectively, 
for all directions ( 0 s). 

We find below that a model spectral function, 
7 rA(k F ,w) = £"/[(w-S ') 2 + £" 2 ], with such a BCS- 
type A (T, 0) indeed reproduces the ARPES spectra, 
whereas the traditionally used assumption of A (T, 0) = 
£peak(F, 0) is invalid. The self-energy (E = E' + iE") 
we use has a minimal representation with two different 
scattering rates, To and Ti [32, 33], 

E(kp, cj) = -iTi + A 2 /[u; + iT 0 ]. ( 1 ) 

Here Ti is a single-particle scattering rate, and it causes 
the broadening of a peak width. On the other hand, T 0 
fills the spectral weight around Fp, and should be viewed 
as the inverse pair lifetime (or pair breaking rate). For 
clarity, we label the former T S i ng i e and the latter T pa i r in 
the rest of this paper. We emphasize that the intensity 
at Ep in a gapped spectrum becomes non-zero only when 
Tpair 7 ^ 0 as simulated in Supplementary Fig. 8 a. Our 
spectra measured at the low temperatures (T <T C ) have 
a negligible intensity at Ep, which ensures that our data 
are almost free from impurity-causing pair breaking ef¬ 
fect. At elevated temperatures, we observe a remarkable 
gap filling (see Fig.Id and Supplementary Figs. 2 and 
5). Significantly, it actually begins from deep below T c , 
which is not expected in a conventional BCS supercon¬ 
ductor. Since the data were measured at the extremely 
high energy resolution (As =1.4 meV), we can rule out 
the possibility, assumed before with setting T pa j r = 0 
[ 6 , 15], that the filling is caused by a spectral broadening 
due to the experimental energy resolution. The inten¬ 
sity at Ep should instead be a signature of intrinsic pair 
breaking, hence it must be taken into account for the 
gap estimation. In passing, we note that the T S i ng i e and 
T pa ir both equally increase the intensity around Ep of 
the momentum-integrated spectrum previously studied 
by ARPES [9, 18, 19] and STM [20, 34] (see Supplemen¬ 
tary Fig. 9), which is therefore incapable of disentangling 
these two different scattering rates. 

Following this consideration, we set T S i ng i e and T pa j r 
to be independent free parameters in Eq.(l). First, we 
performed a spectral fitting to our ARPES data, assum¬ 
ing A (T) = £ pea k(F), which is a traditional way of gap 
estimation (Supplementary Fig. 11a and Fig. 12a). The 
obtained parameter of T sing i e (T) (middle panel of Sup¬ 
plementary Fig. 11a and Fig. 12a) is strongly devi¬ 
ated from a monotonic decrease on cooling, having an 
unrealistic upturn around the temperature at which the 
£peak becomes zero. As already discussed above, this 
anomaly is expected when the spectrum with a single 
peak (£ pea k = 0) has an energy gap (A ^ 0), thus the 
spectral width overestimates the scattering rate. 


We find that this circumstance is corrected by apply¬ 
ing a BCS-type gap function with an onset at 135K and 
89K for OP92K and OD72K (green curves in the bot¬ 
tom panels of Fig.2a and Fig.2b), respectively, regard¬ 
less of the Fermi angle 0. In the Supplementary Fig. 
14 and Fig. 15, we fit Eq.(l) with such a gap function 
A to our ARPES data near the node measured over a 
wide temperature range. As an example, the result at 
0 = 13.5° for OP92K is shown in Fig.5c. The fitting 
curves (red curves) almost perfectly reproduce the data 
(black curves). The obtained T S i ng i e (T)s (Fig. 5a and 
5b) in the gapped region agree with T S i ng i e (T) at the 
node, which can be determined simply from the spec¬ 
tral width. Similarly, the T pa i r (T) curves are also almost 
identical for all the 0s. The consistency in our results 
pointing to the isotropic scattering mechanism validates 
our model spectral function characterized by the BCS- 
type A(T, 0). The famous “hot spots”, at which the 
scattering rate is abruptly enhanced, should be situated 
at much higher 0s (Supplementary Fig. 13). The ap¬ 
plied onset temperatures are almost the same as those of 
Nernst and diamagnetic effects [25, 26], that are viewed 
as signatures of phase-incoherent superconductivity. The 
comparable temperatures are also obtained by the spe¬ 
cific heat measurements [35] and the other spectroscopic 
techniques [36, 37]. Therefore, we assign ~ 135K and 
~ 89K to be the onset temparature of pair formation 
(T p air) of OP92K and OD72K, respectively. This is fur¬ 
ther supported by the signature of pairing seen in the 
behavior of T s i ng i e (T) (Fig.5a and 5b, and Supplemen¬ 
tary Fig. 16); the decrease of its value upon cooling is 
accelerated across T pa i r , showing a deviation from the lin¬ 
ear behavior [38]. The different experimental techniques 
could have different sensitivities to the superconducting 
fluctuation above T c , and actually the terahertz spec¬ 
troscopy estimates a slightly lower temperature scale ( 10 - 
15K above T c ) [39, 40]. Nonetheless, we stress that the 
view that the point-nodal pairing survives above T c is 
compatible in these observations. The doping variation 
of T pair /T c (1.47 and 1.24 for OP92K and OD72K, re¬ 
spectively) obtained in our studies is consistent with the 
phase-fluctuating superconductivity, which merges to the 
superconducting dome with heavily overdoping [41]. The 
competing ordered-phase is, in contrast, claimed to ter¬ 
minate at zero temperature inside the superconducting 
dome [42], thus disagrees with the present gapped states 
observed above T c even in the overdoped sample. 

The relationship between T pair and T sing i e should pro¬ 
vide rich information relevant for the pairing mechanism. 
Intriguingly, the superconductivity occurs when the mag¬ 
nitude of r pa i r is reduced smaller than that of T S i ng i e ; the 
T c is coincident with the temperature at which T S i ng i e (T) 
and r pa i r (T) crosses (magenta circles in Fig. 5a and 5b), 
which provides a simple formula of T sing i e (T c )= r pair (T c ). 
The magnitude of T pair is reported to be comparable 
(120K^150K) among different cuprate families with sig- 
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nificantly different T c s [36]. The r s i ng i e (T) also seems 
to be less sensitive to the different compounds [43, 44]. 
Therefore, the pair breaking effect, which controls the ful¬ 
fillment of T pair < T s i n gi e , is predicted to be a critical fac¬ 
tor determining the T c value of cuprates. Notably, a re¬ 
markable difference in filling behaviors of the spectral gap 
is indeed observed between Bi2212 and Bi2Sr2Cu06+<5 
with about three times different T c s [9, 36]. 

We summarize our conclusion in Fig.6 by drawing a 
schematic temperature evolution of the pairing gap. As 
demonstrated elsewhere [9], the pseudogap competing 
with the pairing disappears far off the node (blue area). 
In this article, we have investigated the nodal region of 
optimally and overdosed samples. We revealed that the 
point-node state persists up to far above T c , against the 
previous expectation [15-17, 20, 21, 24], following the 
BCS-type gap function with the onset at T pair 1.5T C 
in the optimal doping). The ARPRS spectra are repro¬ 
duced by the spectral function with a minimal model, 
which has a single energy gap all the way up to the gap 
closing temperature. It is consistent with the expected 
“pure” pairing state with no contamination by the pseu¬ 
dogap around the node at least in the optimally and over¬ 
dosed regions. While the gap evolution across T pa i r might 
be reminiscent of a phase-transition phenomenon, it is 
against the specific heat data, showing a gradual varia¬ 
tion with temperature above T c [35]. The crossover-like 
behaviors in cuprates should come from the significant 
pair breaking effect, which dramatically enhances across 
T c . To fully understand the present results, insight of 
the spacially inhomogeneous state [45] would be essential. 
However, only that cannot explain our data, since the lo¬ 
cal density of states itself has the behavior of gap filling 
at Ey with temperature [34]. The competing nature of 
pseudogap state evolving around the antinode [1, 3] is 
a plausible source for the unique scattering mechanism, 
which strongly suppress T c from T pa i r . To evaluate this 
speculation, however, the more detailed theoretical in¬ 
puts are required. 

Methods 

Optimally doped Bi 2 Sr 2 CaCii 2 08+<5 (OP92K) and 
overdoped (Bi,Pb)2Sr 2 CaCii208+ ( 5 (OD72K) single crys¬ 
tals with T C =92K and 72K, respectively, were grown by 
the conventional floating-zone (FZ) technique. A sharp 
superconducting transition width of K (OP92K) and 
^3 K (OD72K) were confirmed (see Supplementary Fig. 
!)• 

ARPES data was accumulated using a laboratory- 
based system consisting of a Scienta R4000 electron ana¬ 
lyzer and a 6.994 eV laser. The overall energy resolution 
in the ARPES experiment was set to 1.4 meV for all the 
measurements. In order to accomplish the temperature 
scan of spectra at a high precision, we applied a technique 
of the local sample heating, which thermally isolates the 
sample holder with a heat switch from the lest of the 
system at elevated temperatures. It minimizes the de¬ 


gassing, allowing us to keep the chamber pressure better 
than 2 x 10“ 11 torr during the entire temperature sweep¬ 
ing; no sample aging was confirmed (Supplementary Fig. 
4). This method also prevents the thermal expansion of 
sample manipulator, and it enables us to take data in 
fine temperature steps with automated measurement of 
temperature scan from precisely the same spot on the 
crystal surface, which was essential to achieve the aim of 
the present study. 
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FIG. 1: Temperature evolution of ARPES spectra in the nodal region, a, Dispersion maps at several temperatures 
measured along a momentum cut close to the node in OP92K (a red line in the inset of d). Each map is divided by the Fermi 
function at the measured temperature, b, Temperate evolution of EDCs at kF (a circle in the inset of d) from deep below 
(10K) to much higher than T c (130K). Each spectrum is divided by the Fermi function at the measured temperature, c, The 
same data as in b, but symmetrized about Ep. d, The same data as in c plotted without an offset. The inset represents the 
Fermi surface. The bold orange line indicates the momentum region where the Fermi arc was previously claimed to emerge at 
T c . e, Peak energies of spectra in b and c plotted as a function of temperature, £ P eak- The solid and dashed blue curves show 
the BCS gap function with an onset at T c (92K) and slightly above T c , respectively. Error bars in e represent uncertainty in 
estimating the spectral peak positions. 
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FIG. 2: Disagreement between the ARPES results and the conventional BCS gap function. a,b, Temperature de¬ 
pendence of spectral peak energy, £ pea k(T), at various kp points (color circles in the insets) for OP92K and OD72K, respectively. 
In the bottom of each panel, the same curves of £ pea k (T) are normalized to the maximum value at the lowest temperature. A 
red dashed curve and a green solid curve are the BCS gap function with the onset at T c and T pa i r , respectively (T pa i r =135K 
for OP92K and T pa i r =89K for OD89K). The inset shows the Fermi surface with measured kF points (color circles). The bold 
orange line indicates the momentum region where the Fermi arc was previously claimed to emerge at T c . Error bars in a and 
b represent uncertainty in estimating the spectral peak positions. 
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FIG. 3: Momentum variation of spectra showing the absence of Fermi arc at T c . a, Dispersion maps at T c along 
several momentum cuts (color lines in the inset of d) measured for OP92K. Each map is divided by the Fermi function at 
the measured temperature. The described 0s are the directions of measured kF points (0 is defined in the inset of d). b,c, 
Symmetrized EDCs at kp over a wide range of angle 0 (color circles in the inset of d) at 10K and T c (=92K), respectively. 
The red dotted lines are guide to the eyes for the gap evolution, d, Fermi angle 0 dependence of peak energies of spectra in b 
(^peak(lOK)) and c (£ pe ak(Tc))- The orange arrows indicate the momentum region where the Fermi arc was previously claimed 
to emerge at T c . 
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FIG. 4: Evidence for the point nodal gap surviving above T c revealed by means of the energy shift of spectral 
leading-edge, a, The temperature evolution of ARPES spectra (EDCs) at the gap node for OP92K, normalized to the peak 
intensity of each curve, b-d The same spectra as a, but measured off the node (color circles in the inset of e). e-h Temperature 
dependence of the spectral leading-edge, £le (T), defined as energies where the spectral peaks become half in intensity (marked 
by white arrows in a-d). The bottom panels of f-h plot the difference of £le(T) from the T-linear behaviors (red dashed lines), 
which are fit to £le (T) at high temperatures. The onset temperaturess of the deviation from T-linear behavior are around 135K 
for all the kF points off the node (yellow dashed circles). A kink is seen in the difference curve for f (small black arrow) because 
the thermally populated Bogoliubov peaks in the unoccupied side significantly affect the shapes of leading-edge specially at 
kps close to the node with small gaps. The inset of e shows the Fermi surface with the examined kp points. The bold orange 
arrow indicates the momentum region where the Fermi arc was previously claimed to emerge at T c . 
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FIG. 5: Extracted parameters of the minimal model spectral function required to reproduce the ARPES 
spectra, a, The BCS-type gap function used for the fitting to the OP92K data (a green curve), and the obtained single 
particle scattering rate (r s i ng i e ) and the pair breaking rate (r pa i r ) in Eq.(l). b, The fitting results same as in a, but for the 
OD72K data. The values of r pa i r at high temperatures are not plotted, since the spectral shape is insensitive to the T pa i r 
when A is small or zero, and thus it is impossible to determine the value. A small hump seen in the T pa i r around 75K for 
<j) • 11.7° and 13.5° in a comes from a slight difficulty of fitting to the spectra with a peak-dip-hump shape due to the mode 
coupling, which appears below T c and gets pronounced with approaching the antinode. Magenta circle marks the crossing 
point of T single(T) and T pa i r (T). Large blue arrow indicates the temperature, at which the r s i ng i e (T) deviates from a T-linear 
behavior upon cooling, c, ARPES spectra (black curves) and fitting results (red curves) providing the parameters for 0 = 13.5° 
in a. Error bars in a and b represent uncertainty in fitting the model spectral function to the ARPES spectra. 
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FIG. 6: Schematic pairing-gap evolution based on our ARPES results. a, Temperature dependence of the point 
nodal d -wave pairing gap with BCS-type temperature dependence (inset curve) regardless of directions (0s) along the Fermi 
surface. Temperatures for each curve are indicated in the inset with colored circles. The gapped Fermi surface with a point 
node below T c (b) persists beyond T c (c) up to the temperature of pair formation (T pa i r ). d, Emergence of the gapless Fermi arc 
centered at the node due to the pseudogap evolution around the antinode [9, 27]; while the antinodal region is not observable 
at the low photon energies like 7eV, the studies with higher energy photons have demonstrated that the competing pseudogap 
state emerges at |0| > 25° in the optimal doping [9] . 
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Supplementary Note 1: Samples and Experi¬ 
mental method 

Optimally doped B^S^CaC^Og+s (OP92K) and 
overdoped (B^Pb^S^CaC^Os+tS (OD72K) single crys¬ 
tals with T C =92K and 72K, respectively, were grown by 
the conventional floating-zone (FZ) technique. Magnetic 
susceptibilities for these single crystals are shown in Sup¬ 
plementary Figure 1. Sharp superconducting transitions 
with K (OP92K) and ^3 K (OD72K), indicative of 
a high quality, are confirmed. ARPES data were accu¬ 
mulated using a laboratory-based system consisting of a 
Scienta R4000 electron analyzer and a 6.994 eV laser (the 
6th harmonic of Nd:YV04 quasi-continuous wave with a 
repetition rate of 240 MHz). The overall energy resolu¬ 
tion in the ARPES experiment was set to 1.4 meV for all 
the measurements. 

In order to accomplish the temperature scan of spec¬ 
tra at a high precision, we made a special sample stage; 
we use a heat switch to thermally isolate the sample 
holder with a heater and a temperature sensor (totally 
20 x 30 x 10 mm size) from the rest of the system when 
sweeping the sample temperature. With this technique, 
the samples are efficiently heated up from 10K up to 300K 
while keeping the temperature of manipulator rod other 
than the sample stage lower than 35K. It minimizes the 
degassing and keeps the pressure in the measurement 
chamber better than 2 x 10 -11 torr during the entire 
temperature sweeping. The method of local heating also 
prevents the thermal expansion of the long manipulator 
rod, and thus the sample position becomes unchanged 
with temperature. It enabled us to take data in fine 
temperature steps with the automated measurement of 
temperature scan from precisely the same spot on the 
crystal surface, which were essential to achieve the aim 
of this study. 

In Supplementary Figs. 2 and 3, we show the ARPES 
spectra of OP92K obtained at various kp points with 
sweeping the temperature, with and without an offset, 
respectively. These are the original data used for the 
analysis in the main paper. The ARPES technique 
is sensitive to the condition of sample surface, and 
the biggest challenge in it is to prevent the aging of 


sample surface with time. In Supplementary Figure 
4, we carefully check the stability of sample surface 
in our experiment by comparing two sets of spectra 
measured before and after the temperature scan for 
all the momentum cuts we measured. While it took 
about 12 hours to complete each temperature sweeping, 
almost perfect reproducibility of spectral line shapes 
is confirmed for all the data. This ensures that our 
data measured with the new heating technique are 
reliable and suitable for conducing the detailed analysis 
presented in the paper. Supplementary Figure 5 shows 
the data set similar to those in Supplementary Figure 3, 
but measured for OD72K. 


Supplementary Note 2: Further Evidence for 
absence of the Fermi arc at T c 

In the main paper, we demonstrate that the d -wave 
gap with a point node persists at T c by plotting the 
symmetrized EDCs at various kp points around the 
node (see Fig. 3d). Here we present a further evidence 
that the conclusion is robust and it is not sensitive to the 
selection of kp value. Supplementary Figure 6a shows 
the Fermi-function divided band dispersions (same as 
Fig.3a in the main paper) measured at T c along several 
momentum cuts (color lines in Supplementary Figure 
6c). The corresponding EDCs are plotted in Supple¬ 
mentary Figure 6b. Along the diagonal cut (the most 
left panel), the spectra with one peak crosse Ep due to 
the gap node. In contrast, two peak structure (marked 
with arrows), signifying a gap existence, is observed at 
kp (painted with red) for the other directions. This is a 
direct evidence that the d-wave gap with a point node 
persists, thus the Fermi arc is absent at T c . 


Supplementary Note 3: Temperature evolution 
of the “artificial” Fermi arc 

The Fermi arc [1] is one of the most interested fea¬ 
tures in cuprates. Previous studies have proposed that 
the gapless Fermi arc abruptly emerges at T c , and ex- 
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pands linearly with a further increase of temperature [2- 
4]. The arc length was claimed to be extrapolated to 
zero at T = 0 [3, 4]. The so-called “nodal liquid” behav¬ 
ior has been having a huge influence on the condensed 
matter community. 

In the main paper, we argue that the tracking of 
peak positions of the symmetrized EDCs (£ pea ks) un¬ 
derestimates the gap closing temperature, thus provides 
an “artificial” Fermi arc. Here we demonstrate that 
even such an artificial arc behaves quite differently from 
the previous expectation, when high resolution data 
by a laser ARPES are examined. In Supplementary 
Figure 7a and 7a, we estimate the length of arc with 
£peak=0 as a function of temperature for OP92K and 
OD72K, respectively. In the same panels, the nodal 
liquid behaviors are superimposed. Our results are 
significantly off these lines. As discussed in the main 
paper, the spectral line-shape is determined not only 
by the energy gap (A), but also by the mutual relation 
between the single-particle scattering rate (r s in g ie) and 
the pair breaking rate (T pa i r ). Therefore, it is rather 
reasonable that the growth of arc with £ pea k=0 becomes 
arbitrary. 


Supplementary Note 4: Effect of scattering 
rates, T S ingie and r pa i r , on the spectral line shape 

In the analysis, we used the spectral function with a 
phenomenological self-energy (ref. [5]) given by 


A(k, a;) 


1 


£"(k,a;) 


7T [w — e(k) — E'(k, w)] 2 + S"(k, w) 2 ’ 

E(k, io) = — ir s j ng ie + A 2 /[w + e(k) + irpair]- (SI) 


8b, where the spectra with several T pa i r values are calcu¬ 
lated with fixing T s ingie to zero. Importantly, the energy 
position of spectral peak (£ pea k) dramatically decreases 
with an increase of the T pa i r value as a consequence of 
the spectral filling around Fp, and eventually it becomes 
zero. This simulation indicates that the peak energy un¬ 
derestimates the “real” energy gap (£ pea k < A) when the 
magnitude of T pa j r becomes large. 

Here we should point out that the following self-energy 
has been used to estimate the magnitude of A by a fitting 
to the ARPRS spectra in the many previous reports, 

S(k F ,w) 3 -ir sing ie + A 2 /w. (S2) 

This is the same formula as Eq.(Sl) setting T pa j r = 0, 
and thus it always provides one the spectral intensity of 
zero at Fp as discussed above (see Supplementary Figure 
8a). We emphasize that the spectral weight around Fp 
observed in the ARPES data was previously attributed 
to the spectral broadening coming from the finite energy 
resolution in ARPES (for example, see the supplemen¬ 
tary information of ref. [2]). The extremely high energy 
resolution in our equipment enabled us to reveal that the 
spectral weight at E? is not an artifact due to the exper¬ 
imental resolution, but an intrinsic signature of the pair 
breaking effect. This signifies that the Eq.(Sl) is more 
suitable than Eq.(S2) as the model spectral function of 
cuprates, and that the effect of T pa i r must be taken into 
account to properly estimate the value of A from the 
spectral shape. 

In Supplementary Figure 8c, we simulate a special case 
setting r pa i r = T S i ng i e in Eq.(Sl). The self-energy and the 
corresponding spectral function at kp under this condi¬ 
tion are derived as follow; 


Here A is the energy gap, and e(k) the energy dispersion 
with £(kp) =0. r s i n gi e and T pa i r are the single parti¬ 
cle scattering rate and the pair breaking rate, respec¬ 
tively. In Supplementary Figure 8, we examine the ef¬ 
fect of these scattering rates on the spectral line shape 
at kp with fixing A to 10 meV. To make it realistic, 
all the curves plotted are convoluted with a Gaussian 
that has the width of the experimental energy resolution 
(As = 1.4meV). Nevertheless, note that the A^ value is 
so small that the difference in shape from original curves 
is negligible. Figure 8a shows the T S i ng i e dependence of 
the spectral function when the T pair is fixed to zero. One 
finds that the intensity at Ep is always zero regardless of 
the r s i ng i e value, whereas the spectral width gets broad¬ 
ened with increasing T S i ng i e . This situation strongly con¬ 
trasts to that in our ARPES data (Supplementary Fig¬ 
ure 2 ), which clearly show the filling of spectral weight 
at E p with increasing temperature. The behavior of gap 
filling can be produced by inputting non-zero T pair values 
in Eq.(Sl). It is demonstrated in Supplementary Figure 


E(kp, u j) — — ir T A 2 /[uj T iT], 

A(k F ,w) = X - 

where T = T sing i e = T pair . This formula, consisting 
of double Lorentzian functions, is familiar as the phe¬ 
nomenological BCS spectral function [6-8]. The calcu¬ 
lated spectra for various T values are plotted in Supple¬ 
mentary Figure 8c. One can confirm a couple of discrep¬ 
ancies in these from our data. First, the spectral peak 
in the ARPRS data has an asymmetric shape about the 
peak energy, differently from the Lorentzian shape. Sec¬ 
ondly, the effect of gap filling in this simulation is quite 
sensitive to the magnitude of T, and it is far more drastic 
than that observed in our data at elevated temperatures. 
These lead us to conclude that the Eq.(S3) is not a proper 
formula of spectral function for cuprates either; that is, 
the assumption of T pa j r = T S i ng i e is not suitable. 


+ 


{uj- A) z + r 2 (cj + a)" + t 2 


.(S3) 
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We also stress that the Dynes function [9] given by, 


^dos(<^) = Re 


ijj — iV 

\J{u- ir) 2 - A 2 


(S4) 


is the density-of-states version of Eq.(S3) with Y = 
Tsingie = r pa ir, hence it is also not adequate for cuprates. 
This is understandable because there is no good physi¬ 
cal reason that these two different scattering rates should 
be equal. Our results demonstrate that a proper combi¬ 
nation of T s i n gi e and Tp a i r is required to reproduce the 
data. 

In the previous reports [10-12], the spectrum inte¬ 
grated along a selected momentum cut was utilized to 
study the gap evolution with temperature and its rela¬ 
tion with the electron scattering. Here we argue the ne¬ 
cessity of investigating the one-particle spectrum, rather 
than its momentum integration, for such a study. In Sup¬ 
plementary Figure 9, we demonstrate that a momentum- 
integrated spectrum, Idos(^) = f A(k,uj)dk for A(k, u;) 
of Eq.(Sl), is equally sensitive to the r s i ng i e value (Sup¬ 
plementary Figs. 9a and 9b) and the r pa i r value (Supple¬ 
mentary Figs. 9c and 9d); these two scattering rates both 
fill the spectral intensity around Ep, while the gap size 
A is fixed to a constant value (15 meV). The spectra of 
^dos(^) calculated under the condition of r s i ng i e = r pa i r 
(Supplementary Figs. 9e and 9f) reproduce those in Sup¬ 
plementary Figs. 9a and 9b. As naturally expected, only 
the half values of r sing i e and r pair are sufficient to cause 
the same effect in this case. Notably, the obtained curves 
are identical to Dynes function in Eq.(S4) as pointed 
out above, which is demonstrated for several values of 
T in Supplementary Figure 9g. These simulations pro¬ 
viding identical results (Supplementary Figs. 9a, 9c, 9e, 
and 9g) signify that the momentum-integrated spectrum 
is not capable of extracting the mutual relationship be¬ 
tween the two scattering rates (r s i ng i e and r pa i r ). 

In contrast to it, the effects of r s i ng i e and r pa i r on the 
line shape of the momentum-resolved spectrum A(k, u;) 
are independent, broadening the width and filling the 
weight around Ep, respectively. This circumstance is 
demonstrated in Supplementary Figure 8, and also seen 
in Supplementary Figure 9h, at which the spectra at 
kp are extracted from images in Supplementary Figs. 
9b, 9d, and 9f. (Spectral colors of red, blue, and 
green are corresponding to those of thick dashed lines 
on the images.) Therefore, the precise investigation 
for the line shape of one-particle spectra is required 
to separate the three important parameters (r s i ng i e , 
r pa i r , and A) in Eq.(Sl). The main challenge in this 
task, which has been troubling the ARPES people, is 
to accomplish both of a sufficient energy resolution 
and an ultra-high statistics in the ARPES spectra with 
no noise at a time, which are usually conflicting with 
each other. We have overcome this difficulty by using 


a laser ARPES with a unique cold finger (mentioned 
in ’’Experimental Method”), which was constructed in 
our lab. We emphasize that the Eq.(Sl) is a minimal 
model to reproduce the ARPES data, thus the relation 
among the three physical parameters, revealed in the 
present work, provides a key ingredient to formulate the 
mechanism of the high-T c superconductivity in cuprates. 


Supplementary Note 5: Signature of a gap 
existence in the off-nodal, single-peak spectra 
above T c 

Here we demonstrate that the off-nodal spectra above 
T c exhibit a signature of the gap opening even after be¬ 
coming one-peak structure. The bottom panels of Sup¬ 
plementary Figure 10a and 10b plot the symmetrized 
EDCs of OP92K above T c measured at the node (0 = 0°) 
and off the node (</> = 13.5°), respectively. Two peaks 
seen in the off-nodal spectra merge to a single peak at 
T one =103K (Supplementary Figure 10a). For clarity, the 
spectra above the T one are indicated with rainbow colors, 
whereas those at the lower temperatures with gray colors 
in the both of Supplementary Figure 10a and 10a. 

We find that the single-peak spectra off the node 
(rainbow curves in Supplementary Figure 10a) become 
sharper with increasing temperature up to ~ 135K 
(T p air) before eventually broadened at higher tem¬ 
peratures. The anomalous feature is more clearly 
demonstrated in the upper panel of Supplementary 
Figure 10a, where the spectra below and above the 
T pa i r are separated with an offset. As seen in the the 
magnified insets, the spectral tail shifts to lower energies 
in the range of T < T pair (a magenta arrow), while it 
to the higher energies in T > T pa i r (a black arrow). It 
strongly contrasts to the case at the node (Supplemen¬ 
tary Figure 10b); the spectral tail keeps shifting toward 
higher energies at elevated temperatures in the both 
temperature ranges (black arrows). We summarize these 
results in Supplementary Figure 10c, by plotting the 
spectral width As (HWHM) at and off the node as a 
function of temperature. The anomalous decrease of As 
is seen only for the data off the node (orange circles). As 
discussed in the main paper, a plausible explanation for 
it is that the off-nodal spectra have small energy gaps, 
even though they exhibit single peak shapes, and they 
become sharper up to the temperate (T pair ) at which the 
energy gap completely closes. 


Supplementary Note 6: Spectral fitting under 
the assumption of A (T) = £ pea k (T) 

In the main paper, we claim that the energy of spec¬ 
tral peak underestimates the energy gap (£ pea k < A) 
close to the node at high temperatures, at which the 
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spectral peak width becomes larger than the magnitude 
of the energy gap. Here we perform a spectral fitting 
to our ARPES data with Eq.(Sl) [same as Eq.(l) in 
the main paper], assuming that these two values are 
identical [A (T) = £ pea k (T), upper panel of Supplemen¬ 
tary Figure 11a (OP92K) and Supplementary Figure 
12a (OD72K)], and demonstrate that the extracted pa¬ 
rameter has an unrealistic behavior, reflecting the irrel¬ 
evance of A (T) = £ P eak (T). The middle and bottom 
panels of Supplementary Figure 11a plot the T sing i e (T) 
and rpair(T) of OP92K, respectively, obtained by the fit¬ 
ting to the spectra in the momentum region, where the 
Fermi arc was previously claimed to appear at T c (or¬ 
ange arrow in the insets of Supplementary Figure lib). 
For a comparison, we also plot in Supplementary Fig¬ 
ure lib (OP92K) the corresponding results for A(T) of 
the BCS-like gap function, which are presented in the 
main paper (Fig. 5a). The difference between the two 
fitting-results is very clear; the former (Supplementary 
Figure 11a) shows an abnormal upturn in the T S i ng i e (T) 
curves on cooling in the gapped momentum region, and 
it contrasts to the latter result (Supplementary Figure 
lib) showing almost identical curves in T S i ng i e (T) for all 
0s. In Supplementary Figure 12, we present the same 
data set as in Supplementary Figure 11, but for OD72K. 
Similarly to the case of OP92K, an abnormal upturn is 
obtained in the T S i ng i e (T) when A (T) = £ pea k (T) (the 
middle panel of Supplementary Figure 12a), and it is 
corrected by setting A(T) to be the BCS-like gap func¬ 
tion (the middle panel of Supplementary Figure 12b). 

The behavior of the abnormal upturn is interpreted as 
follows. An energy gap is still open above T c even though 
the spectrum has a single peak, and thus the spectral 
peak width overestimates the scattering rate. Since the 
degree of the overestimation gets higher with an increase 
of the energy gap on cooling, the upturn in the T s i ng i e (T) 
appears. The upturn seems to be more enhanced with 
getting away from the node (see middle panel of Sup¬ 
plementary Figure 11a and Supplementary Figure 12a). 
This is also expected as the overestimation becomes more 
serious toward the antinode with a larger energy gap. 
This scenario is further validated by the fact that the up¬ 
turns start around T ~ 130K and T ~ 90K for OP92K 
and OD72K, respectively, which are consistent with the 
onset temperatures of pair formation determined in the 
main paper. 

One might think that the cuprates are known to have 
a significantly anisotropic scattering mechanism, thus 
the situation in Supplementary Figure 11a (OP92K) and 
Supplementary Figure 12a (OD72K) with the strong 
momentum variation could be intrinsic. To argue against 
this, here we demonstrate that even the optimally doped 
Bi2212, in fact, has an isotropic scattering mechanism 
around the node. Supplementary Figures 13a and 13b 
show the spectra at various kp points over a wide 0 
angle measured at the lowest temperature (T=10K) and 


above the pairing temperature (T=150K), respectively. 
At these temperatures, the value of Tsingie can be esti¬ 
mated simply from the peak width of spectra, since the 
magnitude of T pair is negligible at 10K, and the second 
term in Eq.(Sl) with T pa j r is irrelevant above the pairing 
temperature. The values of T S i ng i e estimated from the 
peak width of spectra in Supplementary Figure 13a and 
13b are plotted in Supplementary Figure 13d with blue 
and red circles, respectively. The magnitude is almost 
constant over a wide 0 centered at the node for both 
the temperatures, which indicates that the scattering 
mechanism near the node is isotropic. This strongly 
supports our assertion that the anisotropic behavior of 
Tsingie with an upturn on cooling is an artifact. Beyond 
|0| ~ 15°, the scattering rate abruptly increases. This 
signifies the evolution of the competing pseudogap, 
which becomes dominant around the antinode [13-16]. 

Supplementary Note 7: Fitting curves repro¬ 
ducing the ARPES spectra 

In the main paper, we demonstrate the spectral fitting 
with a minimal model of Eq.(Sl) to the data of OP92K 
at 0 = 13.5° (Fig.5c), as a typical example. Here we 
present all of the fitting results for the ARPES spectra 
of OP92K and OD72K, which provide the physical pa¬ 
rameters plotted in Fig. 5a (Supplementary Figure 14a) 
and Fig. 5b (Supplementary Figure 15a), respectively. 
Supplementary Figure 14b and 15b shows the ARPES 
spectra (black curves) and fitting results (red curves) for 
OP92K and OD72K, respectively. The data are almost 
perfectly reproduced by the fitting curves. 


Supplementary Note 8: The pairing signature 
detected in the nodal T s i ng i e (T) 

The single-particle scattering rate (T S i ng i e ) at the 
node is extracted simply from the spectral width, thus it 
can be precisely determined without any models. Sup¬ 
plementary Figure 16a and 16b plots the temperature 
dependence of EDC width A^ at the nodal kp point of 
OP92K and OD72K, respectively (same as T s i ng i e (T) at 
0 = 0° in Figs. 5a and 5b). We find that the curves of 
A s{T) deviate from the T-linear behavior (Aeii near (T)) 
in the both samples. It is demonstrated more clearly in 
the bottom panels by plotting the amount of divination 
[A£i inear (T) — A^data(T)]. The onset temperatures of 
the deviation are estimated to be ~ 135K and ~ 89K for 
OP92K and OD72K, respectively (yellow circles). The 
suppression of the scattering rate should be tied to the 
electron pairing. We found that the onset temperatures 
are indeed coincident with the T pair values (gap opening 
temperatures) estimated from the off-nodal spectra in 
the main paper. 
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Supplementary Note 9: Fitting results with a 
different scheme 

The model spectral function of Eq.(Sl) [same as 
Eq.(l) in the main paper] has three parameters of 
the energy gap (A), the single-particle scattering rate 
(Tsingle )5 and the pair breaking rate (r pa i r ). In the main 
paper, we assume the A to be BCS-type gap function 
with the onset temperature much higher than T c , and 
extract the other two parameters of r s i ng i e and r pa i r 
by fitting the model function to the data. Here we 
perform a new analysis with a different fitting scheme 
to justify the BCS-type gap function assumed. In the 
new analysis, we use the nodal r sing i e (T), which can 
be determined without models, for the fitting of the 
off-nodal spectra, and set the other two parameters of 
A and r pa i r free. The underlying idea for it is that the 
scattering mechanism is rather isotropic near the node. 
To make the fitting realistic, we slightly shifted the nodal 
^single(T) in the magnitude. The fitting results (OP92K, 
0 = 11.7°) are demonstrated in Supplementary Figure 
17a. The fitting spectra (red curves) reproduce the data 
(black curves) reasonably well. The extracted A and 
T pa ir also more or less agree to the results presented in 
the main paper, while these values above ~ 120K are 
not properly extracted owing to the small magnitudes 
of A. This result justifies that the A(T) in cuprates has 
the BCS-type gap function with an onset higher than 
T c , and validates the requirement of pair breaking rate 
(T pa ir) to reproduce the ARPES spectra. 
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Supplementary Figure 1: Magnetic susceptibility for single crystals of optimally doped B^S^CaC^Os+s (OP92K) and 
overdoped (Bi,Pb)2Sr2CaCu208+s (OD72K) with an onset T c of 92K and 72K, respectively, used for the ARPES measurements. 
The crystal structure of Bi2212 is drawn in the inset. 
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Supplementary Figure 2: Temperature evolution of ARPES spectra (symmetrized EDCs) at various kps for OP92K (color 
circles in the inset). The corresponding 0 angles (defined in the inset) are described in each panel. The orange arrow in the 
inset indicates the momentum region where the gapless Fermi arc was previously claimed to emerge at T c . 
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Supplementary Figure 3: Temperature evolution of ARPES spectra (symmetrized EDCs) at various kps for OP92K: the 
same data as in Supplementary Figure 2, but plotted with an offset. The top left panel plots the data at the node (0 = 0°). 
Toward the bottom right, the measured kF approaches the antinode. The blue, green, and red curves in each panel are the 
spectra at 10K, T c (=92K), and 150K, respectively. 
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Supplementary Figure 4: Aging check of the sample surface after temperature scan. The upper left panels in a show the 
data along the node (0 — 0°). Toward the bottom right in m, the measured momentum approaches the antinode, n, Fermi 
surface with the measured momentum cuts and kF points. 
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Supplementary Figure 5: Temperature evolution of ARPES spectra (symmetrized EDCs) at various kps for OD72K (color 
circles in the inset). The corresponding 0 angles (defined in the inset) are described in each panel. The orange arrow in the 
inset indicates the momentum region where the gapless Fermi arc was previously claimed to emerge at T c . 
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Supplementary Figure 6: Evidence for the absence of the Fermi arc at T c . a, Fermi-function divided band dispersions at 
T c along several momentum cuts (color lines in c). Same data as Fig.3a in the main paper, b, The corresponding EDCs. The 
spectra at kps (color circles in c) are painted with red. Peak positions are indicated with red arrows, c, The Fermi surface. 
The orange arrow indicates the momentum region where the gapless Fermi arc was previously claimed to emerge at T c . 
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Supplementary Figure 7: The “artificial” Fermi arc. The length of arc where the symmetrized EDCs have a single-peak 
(£ P eak=0) for OP92K (a) and OD72K (b). The nodal liquid behavior, which has a T-linear fashion and is extrapolated to zero 
at T = 0, is superimposed in each panel. 
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k F ,a >) - -iT single + A 2 / [&> + /r pair ], 
A = lOmeV 


® (^single; rpair) — (X, 0) ^ (r s ing^ Tpair) — (0, X) ® (r s ing^ Tp a ir) — (X, X) 



Supplementary Figure 8: Simulation of spectral function in Eq.(Sl) for examining the effect of two different scattering rates 
(rsingle and rpair) on the spectral shape. In all the curves presented, the gap magnitude is set to lOmeV (A = 10 meV). a, 
The Tsingle dependence when r pa ir = OmeV. b, The T pa i r dependence when r s i ng i e = OmeV. c, The spectral shape for various 
values of r s i ng i e and T pa i r , which are set to be equal. In each bottom panel of a, b, and c, all spectra in the upper panels 
are superimposed. While all the curves plotted are convoluted with a gaussian that has the width of the experimental energy 
resolution (As — 1.4meV), the As value is so small that the difference in shape from original curves is negligible. 
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Supplementary Figure 9: Simulations for examining the effect of r s i ng i e and r pa ir on the momentum-integrated spectra. The 
magnitude of energy gap A is fixed to 15 meV for all the cases, a, c, e, Momentum integrated spectra, JdosH = / A(k,uj)dk, 
for the dispersion images in b, d, f, respectively, b, d, f, Dispersion images for one-particle spectral function, A(k,cj), in 
Eq.(Sl) with a linear e(k) (light blue dashed lines). r s i ng i e and r pa ir dependences are studied in (a, b) and (c, d), respectively. 
The case under the condition of r s i ng i e = T pa i r is simulated in (e, f). g, The spectra of Dynes function with several values of 
T. h, Spectra at kF, ^4(1 cf,Cl;), extracted from images of b, d, f. Each color (red, blue, and green) in curves corresponds to 
that of thick dashed line added in the images of b, d, f. 
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Supplementary Figure 10: Signature for the gap existence in the single-peak spectra off the node. Symmetrized EDCs 
of OP92K at T > T c measured off the node (bottom panel of a) and at the node (bottom panel of b). The kp points are 
marked with circles in the inset of a (0 = 13.5°) and b (0 = 0°). Two peaks in the off-nodal spectra merge to one peak at 
103K (Tone)- c, The spectral width (Ae) at T > T one are plotted as a function of temperature. The off-nodal A e(T) decreases 
with temperature and reaches the minimum at ~ 135K (T pa i r ) before it eventually increases at elevated temperatures. The 
anomalous behavior contrasts to the nodal Ae(T), which monotonically increases over the whole range of temperatures. In the 
upper panels of a and b, the spectra above T one with a single-peak (rainbow colors) are extracted. In addition, the spectra 
at T > Tp a i r and T < T pa i r are separated with an offset (black arrow). The insets show the magnified areas marked with the 
dashed squares. The black and magenta dashed arrows indicate the energy shifts of spectral tails with increasing temperature. 
The inset of c plots the energy positions of symmetrized EDCs (£ pea k) as a function of temperature. The green curve indicates 
the BCS-type gap function with the onset of 135K, which fits well the £ pea k(T) at the low temperatures. 
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Supplementary Figure 11: Comparison of fitting results with Eq.(Sl) for OP92K data between two different settings to the 
energy gap, A. a, The A is defined to be equal to the peak position of spectra [A(T) = £ pea k(T)], shown in the top panel, 
b, The A is defined to have a BCS-type gap function with an onset at T pa i r =135K, shown in the top panel. The middle and 
bottom panels plot the obtained fitting parameters, r s i ng i e and T pa i r , at several kF points in the momentum region, where the 
gapless Fermi arc was previously claimed to emerge at T c (see circles and an arrow in the inset of b). The values of T pa i r at 
high temperatures are not plotted, since the spectral shape is not sensitive to the T pa i r when A is small or zero, and thus it is 
impossible to precisely determine the value. 
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Supplementary Figure 12: Comparison of fitting results with Eq.(Sl) for OD72K data between two different settings to 
the energy gap, A. a, The A is defined to be equal to the peak position of spectra [A(T) = £ pea k(T)], shown in the top panel, 
b, The A is defined to have a BCS-type gap function with an onset at T pa i r =89K, shown in the top panel. The middle and 
bottom panels plot the obtained fitting parameters, r s i ng ie and T pa i r , at several kp points in the momentum region, where the 
gapless Fermi arc was previously claimed to emerge at T c (see circles and an arrow in the inset of b). The values of T pa i r at 
high temperatures are not plotted, since the spectral shape is not sensitive to the T pa i r when A is small or zero, and thus it is 
impossible to precisely determine the value. 
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Supplementary Figure 13: Isotropic scattering mechanism around the node, a, b, ARPES spectra of OP92K measured at 
various kF points over a wide 0 angle measured far below T c (T=10K) (a) and above T pa i r (T=150K) (b). In order to remove 
the effect of Fermi cut-off, the curves in a are symmetrized about the Fermi level, and the ones in b are divided by the Fermi 
function at 150K. c, The Fermi surface. The kF points studied in Fig.5a and Supplementary Figure 14 are marked by colored 
circles, d, The peak width of spectra (or r s i ng i e ) in a and b. The orange arrows in c and d indicate the momentum region 
where the gapless Fermi arc was previously claimed to emerge at T c . The almost isotropic Tsingie is obtained in this momentum 
region. 
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Supplementary Figure 14: Fitting results for OP92K data with the model function of Eq.(Sl). a, Same panel as Fig. 5a 
in the main paper: single-particle scattering rate (r s i ng ie) and pair breaking rate (r pa i r ) extracted by setting the energy gap 
(A) to be the BCS-type gap function (a green curve) with the onset of 135K. b, ARPES spectra (black curves) and fitting 
results (red curves) are overlapped. Fermi angle (0) is described in the top of each panel, and the corresponding kF point is 
marked with a circle in the inset of a. We added a small background linear in energy (oc |cj|) to the fitting function A(&f,cj) 
for cf> = 11.7° and 13.5°, in order to properly extract the scattering rates. 
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Supplementary Figure 15: Fitting results for OD72K data with the model function of Eq.(Sl). a, Same panel as Fig. 5b 
in the main paper: single-particle scattering rate (r s i ng i e ) and pair breaking rate (r pa i r ) extracted by setting the energy gap 
(A) to be the BCS-type gap function (a green curve) with the onset of 89K. b, ARPES spectra (black curves) and fitting 
results (red curves) are overlapped. Fermi angle (0) is described in the top of each panel, and the corresponding kp point is 
marked with a circle in the inset of a. We added a small background linear in energy (oc |u;|) to the fitting function A(kF,w) 
for <j) «fi 14.8° and 16.9°, in order to properly extract the scattering rates. 
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Supplementary Figure 16: Signature of the pairing detected in the nodal spectra. Temperature evolution of peak width 
(A^data) of the symmetrized EDCs for OP92K (a) and OD72K (b) are plotted. In the bottom of each panel, the deviation 
of the Asdata(T) from the T-linear behaviors (red dashed lines) are estimated. The onset temperatures of the deviation are 
marked by yellow circles. The T c and T pa i r determined in the main paper are indicated by a light-blue and yellow arrow, 
respectively. 
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Supplementary Figure 17: Fitting results with a different scheme from that used in the main paper; here we input model 
values to r s i ng i e (T) in Eq.(Sl) and extract the other two parameters of A (T) and r pa i r (T). a, The off-nodal ARPES spectra 
of OP92K (black curves) and the fitting results (red curves) at 0 = 11.7° (red circle in the inset of c). We used the nodal 
r s in g ie(T) slightly shifted in the magnitude (orange circles in c), and added a small background linear in energy (oc |u;|) to the 
A(kF,w) in Eq.(Sl) for a realistic fitting to extract the other two parameters (A and r pa ir)- b, The obtained A(T) are plotted 
(light blue circles). For comparison, the BCS-type gap function with the onset of 135K (green curve) and the peak energies of 
symmetrized EDCs (£ pea kS in Fig.2a of the main paper) are also shown, c, The used r s i ng i e (T) for the fitting analysis (orange 
circles) and the extracted r pa i r (T) (green circles) are plotted. For comparison, the r s i ng i e (T) (open black squares) and r pa i r (T) 
(open black triangles) obtained in the main paper with a different fitting scheme (Fig.5a of the main paper) are also shown. 













































































